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Abstract 

We discuss the application of the method of characteristics to the 
renormalization-group equation for the perturbative QCD series within 
the electron-positron annihilation cross-section. We demonstrate how one 
such renormalization-group improvement of this series is equivalent to a 
closed-form summation of the first four towers of renormalization-group 
accessible logarithms to all orders of perturbation theory. 
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The renormalization group (RG) has played a central role in our understand- 
ing of quantum field theory [1-10] especially since the discovery of asymptotic 
freedom [11-14]. The central idea of the renormalization group is the insensi- 
tivity of physical quantities to the mass scale /z 2 introduced in the course of 
regularizing and eliminating infinities within perturbative calculations. Explicit 
dependence of a perturbative series on it 2 is compensated by fi 2 dependence in 
masses and coupling constants characterising that series. Indeed, the replace- 
ment of such quantities by running quantities that are explicitly functions of 
li 2 constitutes what is generally denoted by "RG-improvement" of a pertur- 
bative expression [15]. The numerical value of a calculation to a given order 
of perturbation theory still depends upon the numerical value of /i 2 , entailing 
the introduction of either prescriptions (e.g., mi,/2 < fi < 2nib for semileptonic 
6-decays) or procedures [16, 17] to obtain optimal values of /z 2 . 

However, such substitutions do not in themselves take full advantage of all 
information accessible from the renormalization-group equation (RGE), which 
also determines portions of the perturbative series beyond the order of perturba- 
tion theory to which calculations have been explicitly performed. Application of 
the RGE to one-loop expressions has long been known to determine the leading 
logarithm contribution to each subsequent order of perturbation theory. The 
RGE can similarly be used in conjunction with two-loop calculations to deter- 
mine next-to-leading logarithm contributions to all subsequent orders in pertur- 
bation theory - indeed the application of the RGE to an n th loop perturbative 
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expression is sufficient to determine the contribution of n successively-subleading 
logarithms to all orders in the perturbative expansion parameter. 1 Such RGE 
methods for obtaining and summing "RG-accessible" logarithms to all orders 
of perturbation theory have been applied to effective potentials [19] and actions 
[20, 21], QCD correlation functions [22, 23], QCD contributions to decay rates 
[22], and even the high-energy behaviour of the WW — > ZZ cross-section [22], 
a process dominated by Higgs boson exchanges. A related RG-summation of 
dimensionality poles in the expansion of the bare coupling constant in terms of 
its renormalized analog has been developed in ref. [24] and (for thermal field 
theory) in ref. [25]. 

The point we wish to emphasize is that the summation of higher order log- 
arithmic contributions is quite distinct (and a substantial improvement over) 
what is usually understood to be RG-improvement, the incorporation of run- 
ning masses and coupling constants into perturbative expressions taken to a 
given order. Indeed, such inclusion of all RG-accessible logarithms within per- 
turbative series is forcefully advocated in ref. [26]. Series which incorporate 
summation of RG-accessible logarithmic contributions to all orders of perturba- 
tion theory have been seen to exhibit much less dependence on fj, 2 than series 
which utilize running masses and coupling constants to a fixed calculational 
order [22, 23]. This latter approach, however, devolves from the method of 
characteristics [27], a standard approach to first-order partial differential equa- 

lr The possibility that infrared effects might alter this is addressed in ref. [18] 
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tions such as the RGE [28]. We demonstrate below how this same method 
of characteristics can be extended to obtain summations of leading and three 
successively-subleading towers of logarithms to all orders of the perturbative 
QCD series for the electron-positron annihilation cross-section. 

The total cross-section for e+e~-annihilation, R(s) = <r(e + e~ — ► hadrons)/ 
er(e + e~ — > can be extracted from the imaginary part of the QCD vector- 

current correlation function [29], a perturbative expression that necessarily de- 
pends upon a renormalization mass scale /i: 

R(s) = 3j2Q 2 fS[x(n 2 ),log(n 2 /s)]. (1) 
/ 

The expansion parameter x(p 2 ) = a s (n 2 )/ir is proportional to the running 
QCD coupling constant, and renormalization mass scale \x is a by-product of 
the regularization procedure for indentifying and excising infinities from the 
underlying correlation function, as discussed above. Since R(s) cannot depend 
on this unphysical scale parameter, it follows that 

"'ff 1 — 

where 

/?(z(M 2 ))=M 2 ^ (3) 

with an appropriately chosen boundary condition for (3) [e.g. x(A 2 ) — oo or 
x(M 2 ) = 0.118/7r]. Perturbative series expansions of S[x(p, 2 ), log (^^l an d 
(3(x) are seen to take the form 
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S[x, L] = 1 + ^ ^ T n , m x n L m (4a) 

n—1 ??i—0 

oo 

/?(x) ^-x 2 ^/?^ (4b) 

fe=0 

where a; = x(/i 2 ) an d = l°g(M 2 / s )- Generally /3(x) is determined by relating 
the bare coupling to the renormalized coupling [7, 9, 24], although it can also be 
extracted directly from eq. (2) [20, 21, 30]. Indeed, explicit Feynman diagram- 
matic calculations to four-loop order have determined 0o, (3i, (h, (h [13, 14, 31] 
as well as 7i ; o, T24, T2,o, ?3,2, ^3,1 and T^fi [29], and these results (as tabulated 
in Table I of rcf. [23]) are manifestly consistent with eq. (2). 

However, it is possible to utilize eq. (2) to extract higher-order coefficients 
within S[x,L] than those tabulated in Table 1. It is easily seen [23] that Ti ; o 
and Po determine all leading logarithm coefficients T"„ )n _i for n > 1; similarly 
additional knowledge of T 2 .o and f3\ is sufficient to determine all next-to-leading 
logarithm coefficients T„. n _ 2 f° r n > 2; knowledge of T 3i o and /3 2 permits deter- 
mination of T„,„_ 3 for n > 3, and so forth. In ref. [23], the double summation 
in eq. (4a) is reorganised into the form 

00 

S[x,L] = l + ^x"S„(xL) (5) 

n=l 

where the functions 

00 

S n ( U ) =Y, T n+k,kU k (6) 
k=0 

are completely determined by knowledge of the "RG-accessible" coefficients 
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T n +k,k- One can show that eq. (2) gives rise to a nested set of first order 
differential equations for the functions S n (u): 

^ j^s k = i^Ml £ A + fc - 5„(0) = t„. . 

(7) 

These equations are derived and sequentially solved in ref. [23]. When one 
applies this "RG-summation" to the series (5) within R(s), the dependence of 
R(s) on fj 2 is substantially reduced [22, 23]. This is not surprising, as the exact 
result for R(s) is necessarily independent of ^i 2 [the RGE is just a statement 
of this independence] , and the inclusion of higher-order logarithm contributions 
to R(s) via (5) is expected to approximate the exact result more closely than 
truncation of eq. (4) to a given order. 

As discussed above, the method of characteristics [27,28] provides a comple- 
mentary procedure for obtaining information from the renormalization group 
equation (2). To illustrate this method, consider the first-order partial differen- 
tial equation 

A(x,y)=0 (8) 



where / and g are given functions, and where A(x, y) may be indentified as some 
field-theoretical amplitude characterised by quantities (e.g, coupling constants) 
x and y. If A (x,y) is a solution to eq. (8), then so is A (x(t),y(t)), provided 
that 

§ = /(*(*),»(*)), f t =g(x(t),y(t)) (9) 
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with initial conditions x(0) = x, y(0) = y. One then sees from eqs. (9) that 







A (x,y) = f t A (x(t),y(t)) (10) 



The initial conditions ensure that A (x(t), y{t)) is a solution of eq. (8) when t = 
0. Since eq. (10) implies that A (x(t),y(t)) is independent of t, A (x(t),y(t)) 
is necessarily a solution to eq. (8) for all values of t. Eqs. (9) and (10) provide 
the justification for replacing the variables x and y with their corresponding 
characteristic functions x,y in the amplitude A(x,y). 

For the RGE (2), as applied to the field theoretical series S[x,\og(fi 2 /s)], the 
role of / and g as dependent variables is assumed by /z 2 and /?, in which case 
correspondence to eqs. (9) requires running values for these functions 



dt 
dx(t) 



= n\t), (ii) 

P(x(t))- (12) 



dt 

The usual prescription for RG improvement is to identify t with log(^i 2 ) [i.e., 
ft 2 — e* = n 2 via eq. (11)], in which case eq. (12) becomes eq. (3). Indeed, this 
construction provides the justification for having the coupling constant x run 
with [i 2 within the perturbative series (4a) [27]. 

However, it is entirely valid to let t be arbitrarily chosen in eqs. (11) and 
(12), up to initial conditions fl 2 (0) = ^i 2 , x(0) = x(ji 2 ) that establish contact 
with a known solution to eq. (2). Thus the "running coupling" x(p?) may be 
employed to serve as an initial condition for the characteristic function x(t). 2 

2 The dimensional regularization equation relating the bare (gs) and renormalized (g) cou- 



In order to keep track of the order of perturbation theory to which we are 
working, we follow the approach of ref. [32] by introducing an expansion pa- 
rameter h such that 



t — >t/Ti (13a) 
x — ► xh (13b) 
so that the characteristic equations (9) become 

n 2 ^ = -x 2 h 2 f2x n h n p n (15) 

(XX 

n— 

and the series expansion (4a) becomes 

OO 71—1 

5[x,Z] = l + ^^T„, m ?i"x"Z m (16) 

n=l m— 

[L = log (/Z 2 (i)/s)]. From eq. (14) and the /2(0) = /z initial condition, we see 
that 

jl 2 (t)=fj l 2 e t / h . (17) 



We now express x(t) as a perturbative expansion 

OO 

s(t) = (18) 



n=0 



pling constants is = /*' X]™ o a )i,tf 2fc+1 / f ' [7, 9]. Since gs is a bare parameter 

independent of fi, the renormalizcd coupling-constant g is necessarily a /j-dcpcndcnt quantity, 
i.e., a function of //,. 
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with xo(0) = x(p?) and x n (0) = for n > 0. Upon subsituting eq. (18) into 
cq. (15), we obtain a nested set of linear first-order differential equations for 
the variables x n (t) when n > 0: 



dxi 
~~dT 



^2. = -foxl, x (0) = x{ji 2 ), 



+ (2p x (t))x 1 = -Pi4(t), Si(0) = 0, 



(19) 
(20) 



dx 2 
~~dt 



+ {2p x {t))x 2 = -p 2 xt{t)-3p 1 xl{t)x 1 {t)-p xi{t), x 2 (0)=0, (21) 



dx 3 
~~dt 



+ (20ox o (t)) x 3 



3/3i (x 2 (t)x 2 (t) + x (t)x 2 (t)) 
2(3 x 1 (t)x 2 (t), x 3 (0) = 0. 



(22) 



The solution to eq. (19) is 



x (t) 



x(n 2 



l + (i x{p?)f 

If we substitute eq. (23) into eq. (20), we find the solution to cq. (20) to be 

/?i zV)log(l + pWM 2 )*) 



(23) 



xi(t) 



(24) 



0o (l + A^O 2 )*) 2 
Similarly, subsitution of eq. (23) and (24) into eq. (21) leads to a solution for 

x 2 (t) 



x 2 {t) 



x 3 (fi 2 ) 


\(ft 


&\ 


(l+fox^W 


.U 





ft 
fto 



ft 

fto 



+ ^\og 2 (l + p x^ 2 )t) 



(25) 
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and substitution of eqs. (23), (24) and (25) into eq. (22) leads to a solution of 



(Eqs. (23), (24) and (25) also appear in rcfs. [21,32].) Eqs. (23) - (26) provide an 
expansion (18) of the solution to (15) that is distinct from the usual perturbative 
expansion. 

If we substitute the series (18) into the expansion (16) for S[x,L] we find 
the following solution to the renormalization group equation (2): 



S[x,L] = l + h[Ti t0 x (t)] 

+ h 2 [Ti, ii(t) + (T 2 , + T 2A L)x 2 (t)]] 

+ h 3 [T lfi x 2 (t) + (T 2 , + T 2 ,iL) (2x (t)x 1 (t)) 
+ (T 3 ,o+T 3il L + T 3 , 2 L 2 ) x 3 (t)] 

+ h 4 [Ti, i 3 (t) + (72,o + T 2>1 L){x\(t) + 2x {t)x 2 (t)) 
+ (73,o + T 3 ,iL + T 3 . 2 L 2 ) (3x 2 (t) Xl (t)) 
+ (74,o + T 4A L + T^ 2 L 2 + T 4 , 3 £ 3 ) x 4 (t)] 

+ ... (27) 

Now, if t = 0, we find from eq. (17) that jl 2 = fj, 2 , L — \og([i 2 /s) — L. Since 
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+ ^log 2 (l + /W)0 
z Po 



log 3 (1 + I3 x(p 2 )t) 



(26) 



xo(0) — x{jj, 2 ) and alfc(O) = for fc > 1, we see that eq. (27) recovers the original 
expansion (4a) when K = l. 

It t is a non-zero constant, the solution (27) provides a means for obtaining 
all coefficients T n m with m i ^ in terms of coefficients T^ q. To see this, let 
t = hlogk, in which case we see from eq. (17) that p, 2 — k^, 2 , L = log(k[i 2 /s) = 
L + logfc. If we substitute t = hlogk into eqs. (23) - (26) and note that 

log (1 + h/3 x(iJ, 2 )t) = h(3 x(p 2 )\ogk 

h 2 



y/3 2 xV)log 2 fc 
h 3 



+ y/3 3 xV)log 3 fc + ..., (28) 
we find upon further substitution into eq. (27) that 

S[x,L]\ t =hiogk = l + hx(p 2 )T lfi 

+ h 2 x 2 (p, 2 ) [(T 2i0 -PoTi, logfc) + T 2il log(fc M 2 / S )] 

+ h 3 x 3 (p 2 ) { [T 3i0 - (2T 2 , /3 + T li0 /9i) logfc + T lj0 $ log 2 fc] 

+ [T 3 ,i - 2A)T 2 ,i logfc] log(fc M 2 /s) + T 3 , 2 log 2 (fc M 2 / s )} 
+ & 4 x 4 ( M 2 ) {[T 4 , - (3/? T 3 ,o + 2/?iT 2 , + /3 2 T lj0 ) logfc 
+ + ^o/3iTi, ) log 2 fc - flj}Ti l0 log 3 fc 

+ [T 4 ,i - (3A)T 3 ,i + 2/3^2,0 logfc + 3/3 2 T 2>1 log 2 fc] log(fc^ 2 / s ) 
+ [T 4 , 2 - 3/3 r 3 , 2 logfc] log 2 (fe// 2 /a) +T 4 , 3 log 3 (fc M 2 / s )} 
+0 (ftV(/i 2 )) . (29) 

Now if we rewrite the original series expansion (4a) with L = log(/x 2 /s) = 
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\og(k[i 2 /s) — log A;, we obtain 

S[x,L] = 1 + .t( M 2 )T!,o 

+ x 2 {p 2 ) [(T 2 , -T 2a log fc)+ 72,i log(fcM 2 A0] 
+ x 3 (p 2 ) [(T 3l0 - T 3 ,i logfc + T 3 , 2 log 2 fc) 

+ (T 3 , 1 -2T 3 , 2 logfc)log(fc M 7 S ) 

+T 3 , 2 log 2 (fc M 2 / S )] 
+ x^p 2 ) [(T 4 , - T 4 ,i logfc + T 4 , 2 log 2 k - T 4 , 3 log 3 fc) 

+(T 4 ,i - 2T 4 , 2 log k + 3T 4 , 3 log 2 fc) log(fc M 2 /s) 

+ (T 4 , 2 - 3T 4 , 3 logfc) log 2 (fc M 2 fa) +T 4 , 3 log 3 (fc M 7*)] 
+ 0(*V))- (30) 

When 7i = 1, eqs. (29) and (30) must be equal, since the solution (27) to 
the renormalization-group equation (2) has been constructed... 

1) ...to coincide with eq. (30) at t = 0, 

2) ...to be independent of the choice for t via the method of characteristics 
[eq. (10)], and 

3) ...to be given by eq. (29) when t — % logfc. 

Direct comparison of eqs. (29) and (30) when % = 1 shows that 

72,i = PoTifl, 

T 3 ,i = 2T 2 ,oA>+/?iTi i0 , 



^3,2 


= /3o^i,0) 




24,1 


= 3^)13,0 + 2^2,0 + ^1,0, 




?4,2 


= 3^T 2 ,o + ^o/3iri,o, 




T 4 ,3 


= Po^i.o, 


(31) 



relations that can also be obtained [22] by direct substitution of the series (4a) 
into the renormalization-group equation (2). Thus, the method of characteristics 
is seen to determine all logarithmic coefficients to the order of perturbation 
theory considered. 

However, a more powerful application of the solution (27) occurs by setting 
t = h\og(s/ /i 2 ), ensuring via eq. (17) that jl 2 = s, that L = 0, and that factors 
of 1 + h(3ox([i 2 )t in cqs. (23) - (26) become 1 - (3 x(p, 2 ) log(/z 2 /s) = w in the 
h — ► 1 limit. In this limit, eq. (27) generates the following series: 

S[x,L] = 
+ 
+ 
+ 
+ 
+ 

This series explicitly reproduces the series (5) obtained via the successive solu- 
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l+aV^+^V) 
w 



'-2,0 



*V) 
Pi 

pi ' 

5/3? 



3,0 



2T 2 , -r!- logw 

Po 



PL 

.With 

PI 
W1P2 



Pi 

'Po 

PlP2 
PI 

wlogw 



Pi rp , 

Po 

A 2 . ft 2 ^/3 2 , 

^2 ~ — ) ( w ~ : ) ~ ^2 lo S w + 152 lo ? 



2/?o 



/?o 2 



2T 2,0-32- 



Pl 
'Pi 



PL 

2/3 3 

3T 3 , | 
Po 



77^+3T 2 , ^ log'w 



2po 3 
Po 



w 



+ 



-3 



Po 

Pf 
PI' 

A 

2/?o 
logw 



/3l/?2 

2T 2 , 



Po 
2T2.0 



/5o 2 
# _ P 
P\ 



PJ 
PI 



P2 

Po 



Po 



(32) 



tions to the differential equations (7). The coefficient functions Si(xL), S 2 {xL), 
Ss(xL) and Sa{xL), as calculated in ref. [23], arc reproduced in cq. (32), 
demonstrating how the method-of-characteristics approach to the renormaliza- 
tion group equation can recover the results obtained via summation of leading, 
next-to-lcading, and successively subleading logarithm factors to all orders of 
perturbation theory. Such results as eq. (32) represent the optimal possible RG- 
improvement of the original perturbative series (4a), insofar as they incorporate 
all RG-accessible coefficients of logarithms occurring within that series. 

As a final note, the ambiguity in the choice of k such that t — Ti log k [thereby 
leading to the series of eq. (29)] is equivalent to the ambiguity noted in ref. 
[25]. In [25], this ambiguity was viewed as a consequence of shifting the initial 
condition of eq. (7) to the equation (6) defining S n (u); by replacing log(/z 2 /s) 
with log(fc^ 2 /s) — log(fc) in the series (4a), that series becomes 

oo n— 1 

S[x,L] = l + J2J2 T n,m^ n (LT (33) 

71=1 771=0 

where V = \og(k/j, 2 /s) and where T[ = T 1>0 = 1, T' 2 X = T 2 ,i,T^ = T 2 , - 
T 2 ,i log(fc), etc. The initial condition in eq. (7) is now replaced by S n (0) — T' nS) . 
Of course, when one sums to all orders in perturbation theory, the dependence 
on k within eq. (29) will drop out. 
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